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Abstract 


Computational phenomena (i.e., spurious supersaturation and negative mixing ratio of 
cloud water) usually exist in cloud-resolving models when the time step for explicit 
integration is larger than a microphysical timescale in clouds. In this paper, the 
microphysical timescales in clouds are studied, showing that the timescale of water vapor 
condensation (or cloud water evaporation) is smaller than 10 s — the order of a typical 
time step for cloud-resolving models. To avoid spurious computational phenomena in 
cloud-resolving modeling, it is suggested that moist entropy be used as a prognostic 
thermodynamic variable, and temperature be diagnosed from that and other prognostic 
variables. A simple numerical model with moist entropy as a prognostic variable, for 
example, is presented to show that spurious computational phenomena are removed when 
moist entropy is used as a prognostic variable. 



1. Introduction 


With the continuous increase in computational power, cloud-resolving models will be 
used in the near future to simulate explicitly clouds and large-scale circulations for their 
interaction. For the correct simulation of that interaction, the models should well 
represent not only the processes in relation to the life cycle of individual clouds systems 
but also the processes in relation to large-scale circulations. 

In the Tropics, deep cumulus clouds act as an engine for the atmosphere, driving 
large-scale vertical circulations through convective heating (Riehl and Malkus 1958). 
However, many questions on convective heating need to be addressed through 
observations and modeling (e.g., Simpson et al. 1988, Tao and Adler 2003a). One of the 
questions is on the “efficiency” of the engine. Since clouds consist of small particles, the 
“efficiency” is related to the microphysics in clouds. Hence the microphysics should be 
represented properly in a cloud-resolving model in order to accurately simulate the 
interaction between clouds and large-scale circulations. 

On the other hand, small particles in clouds absorb and emit radiation, which can 
change the atmospheric energy budget and in turn large-scale circulations. As a result, 
they are coupled with large-scale circulations (e.g., Raymond 2000, Raymond and Zeng 
2000). Since large-scale vertical circulations are sensitive to the atmospheric radiative 
cooling rate (e.g., Zeng et al. 2004) and the radiative cooling rate is modulated by clouds 
through cloud microphysics (e.g., Albrecht and Cox 1975; Baker 1997), cloud 
microphysics must be represented properly in a cloud-resolving model for the accurate 
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simulation of radiation and its role in the interaction between clouds and large-scale 
circulations. 

In current cloud-resolving models, there are computational phenomena such as 
spurious supersaturation near cloud edges (e.g., Grabowski 1989), and special adjustment 
techniques are usually introduced to reduce their negative effects, benefiting the 
numerical simulation of individual cloud systems (e.g., Tao et al. 1989). Some 
computational phenomena originate in the choice of prognostic thermodynamic variables 
(Grabowski 1989; Ooyama 1990; Zeng 2001; Zeng, Tao and Simpson 2004). In a new 
trend to avoid them, moist entropy is used as a prognostic variable in modeling, and 
temperature is diagnosed from that and other prognostic variables (Raymond and Blyth 
1986; Ooyama 1990, 2001; Zeng 2001, Zeng, Tao and Simpson 2004). Raymond and 
Blyth (1986) used moist entropy and the total mixing ratio of airborne water (water vapor 
and cloud water) as prognostic variables in a parcel model. Ooyama (1990) analyzed the 
thermal relations of moist entropy and suggested moist entropy be used as a prognostic 
variable in multi-dimensional numerical models. Ooyama (2001) and Zeng (2001) 
constructed two- and three-dimensional models with moist entropy as a prognostic 
variable, respectively, to simulate warm clouds. In their models, neither spurious 
supersaturation nor negative mixing ratio of cloud water appears. Thus it is logical to 
extend moist entropy as a prognostic variable to simulate cold clouds. Zeng, Tao and 
Simpson (2004) derived a precise equation for moist entropy, providing a theoretical 
basis for moist entropy to be used as a prognostic variable in the modeling of cold clouds. 
In this paper, another aspect of moist entropy is addressed: the benefit of using moist 
entropy as a prognostic variable. 
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Spurious computational phenomena, as shown in section 2, exist in a numerical model 
when the time step for explicit integration is larger than a physical timescale. The 
phenomena manifest themselves in cloud-resolving modeling as spurious supers aturati on 
and negative mixing ratio of cloud water. Since the proper simulation of the mixing ratio 
of cloud water is vital in long-term cloud-resolving modeling (e.g., Robe and Emanuel 
1996), the origin of computational phenomena should be explored so that they can be 
removed. Thus, the microphysical timescales in clouds should be studied first in contrast 
to the time step for cloud-resolving modeling. 

The microphysical timescale of water vapor condensation (or cloud water 
evaporation), as shown in section 3, is smaller than 10 s — the order of a typical time 
step for cloud-resolving models. Thus, the small timescale can bring about computational 
phenomena in cloud-resolving modeling. To avoid spurious phenomena, it is suggested in 
section 4 that moist entropy be used as a prognostic variable in the place of temperature. 
That strategy is different from other technologies on cloud condensation modeling (e.g., 
Grabowski 1989; Tao et al. 1989; Grabowski and Smolarkiewicz 1990; Margolin, 
Reisner and Smolarkiewicz 1997). 

This paper is composed of five sections. In section 2, a spurious computational 
phenomenon is illustrated with the aid of a simple numerical model. In section 3, the 
expressions for microphysical timescales in clouds are derived, and the timescales are 
compared with the time steps used in cloud-resolving modeling. In section 4, the choice 
of prognostic thermodynamic variables in cloud-resolving modeling is discussed, and a 
numerical model with moist entropy as a prognostic variable is presented, as an example, 
to show that computational phenomena are removed when moist entropy is used as a 
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prognostic variable. In section 5, a summary is given. Except for special illustrations, the 
paper follows the symbol definitions in Appendix A. 

2. Spurious computational phenomena 

In this section, an example of spurious computational phenomena is illustrated with a 
simple numerical model whose corresponding differential equation is 

d<j>ldt = -0/r , (2.1) 

where the variable ^is a function of time t and the constant ris a timescale. The variable 
(jf= 1 when t=0. Equation (2.1) is solved analytically with <j> = exp (—tit) . 

For explicit integration, Equation (2.1) is discretized as 

-fAt/r (2.2) 

where superscripts indicate time level and A t is the time step for integration. Equation 
(2.2) is solved with <j) n =(1-A f/r)" . Obviously the numerical solution of Equation (2.2) 
is close to the analytical solution of Equation (2.1) when 

A t < r . (2.3) 

Otherwise, the value of / blows up or oscillates spuriously around zero. 

The computational phenomenon when A t>r is shown in Figure 1, where At=10 s 
and the timescale ris 2, 5 and 12.5 s, respectively. The numerical solutions of Equation 
(2.2) and their corresponding analytical solutions are shown with thick and thin lines, 
respectively, in the same pane for comparison. The values for 0 blows up when t=2 s 
(or At > 2r ). When At = 2r = 10 s, <[> oscillates between -1 and 1 spuriously. When 
r=12.5 s (or At < r ), the numerical solution is close to the analytical solution. 
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Many microphysical processes in clouds are governed by Equation (2.1) or similar 
equations with different timescales (see section 3 for details). When a microphysical 
timescale is so small that the condition (2.3) is violated, spurious behavior appears. In the 
next section, microphysical timescales in clouds are studied so as to provide theoretical 
evidences for the construction of cloud-resolving models. 


3. Microphysical timescales in clouds 

a. Condensation growth 

An air parcel with uniform cloud droplets is studied for the timescale of water vapor 
condensation. Only water vapor condensation on cloud droplets or the evaporation of 
cloud droplets exists in the parcel. The concentration and the radius of droplets in the 
parcel are denoted with N c and r c , respectively, and the mixing ratio of cloud water with 
q c . Thus 



The growth rate of a droplet due to water vapor condensation is expressed as (e.g., 
Pruppacher and Klett 1997) 

dr c _ Kiqjq 

vyw — , (3-2) 

dt r c 

where q v is the mixing ratio of water vapor, q vsw the saturation mixing ratio of water 
vapor over water and 

K = <J>£ ' KKT 1 + P.RJ!E„DJ ' . (3.3) 

Differentiating Equation (3.1) with respect to time, and then substituting Equation (3.2) 
into the resulting equation yields 


6 



(3.4) 


~~ = -~T~ = ~^^wAv r c~ S ~(~ 1 ) • 

dt dt p q vsw 

Assume that the parcel is adiabatic and stationary. Thus, the energy equation is written 
approximately as 


r dT - r 

p dt v dt 


(3.5) 


Substituting Equation (3.4) into (3.5) yields 




dt 


P R V C P T< q 


(3.6) 


with the aid of the Clausius-Clapeyron equation 


dhE„ _ L v 
dT RT 


2 ‘ 


(3.7) 


The analogy between Equations (3.6) and (2.1) shows the timescale of water vapor 
condensation 


T = (1 + 


L 2 

vsw ^v \~1 

rr? / 


pq v 


KC p T 


47ip w A w N r r r 


(3-8) 


which is consistent with that of Squires (1952) and Korolev and Mazin (2003). 

Equation (3.8) shows that the timescale is a function of temperature and pressure. An 
ideal case is discussed here to show the variation of the timescale with height. Assume 
that cloudy air has the same temperature and pressure as its ambient air. The ambient air 
is stationary. Its surface pressure is 1013.25 hpa, and temperature decreases linearly with 
height from 288 at z=0 to 216.5 K at z=ll km. For this case, the timescale of water vapor 
condensation against height is computed with Equation (3.8), and the results for N c r c = 500 
and 3000 pm-cm' 3 are displayed in Figure 2 with thin and thick lines, respectively. If 
r c =10 pm, N c r c = 500 and 3000 pm-cm' 3 are representative of marine and continental air, 
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respectively (e.g., Pruppacher and Klett 1997). Both timescales for N c r c = 500 and 3000 
pmcm' 3 , as shown in the figure, are smaller than 10 s, implying that water vapor 
condensation on to cloud droplets (or the evaporation of cloud droplets) is not suitable for 
explicit simulation in cloud-resolving modeling. 


b. Other diffusion growth 

Similar to the derivation for Equation (3.8), the timescale of water vapor deposition on 
small ice particles is obtained as 


r = ( 1+ mv 


R V C P T wan* 


(3.9) 


where q vsi is the saturation mixing ratio of water vapor over ice, r, the radius of small ice 
particles, N, the concentration of the ice particles and 

A = <J>£ / K a RJ 2 + P'Rjr / E si D v r • (3-10) 

Following the procedure in Figure 2, the timescale for vapor deposition is computed 
with Equation (3.9), and the results for A,r,=10' J , 10' 1 and 10 1 pm-cm° are displayed in 
Figure 3. If r,=10 pm, jV,r,= 10' 3 , 10' 1 and 10 1 pm-cm' 3 represent the air with low, 
moderate and high concentrations of ice particles, respectively (e.g., Pruppacher and 
Klett 1997). As shown in the figure, the timescale decreases with the increase of Ni for a 
given r,. However, even when the concentration of ice particles is high (or Nj= 1 000 liter 
’), the timescale is still larger than 10 s, implying that water vapor deposition to small ice 
particles is suitable for explicit simulation in cloud-resolving modeling. 

The timescale of rainwater evaporation can be obtained in steps analogous to Equation 
(3.8). Introducing the ventilation coefficients f v for water vapor diffusion and A for heat 
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transport in Equation (3.2) for raindrop evaporation and assuming f v =fh (e.g., Pruppacher 
and Klett 1997), the time scale of rainwater evaporation is obtained 


T = (1 + 


Ivsw L 1 

KC p t 7 


P9v 


4np w \f v N r r • 


(3.11) 


where N r and r r are the concentration and the radius of rain drops, respectively. Let 
r^=10 3 pm, iVr= 10 3 m' 3 , / v =10, T - 280 K and p=850 hpa. The timescale of rainwater 
evaporation is 323 s. Since the timescale is much larger than 10 s, the order of time step 
in cloud-resolving modeling, rainwater evaporation is suitable for explicit simulation in 
cloud-resolving modeling. 

Similarly, the time scale of water vapor deposition on precipitating ice particles is 
obtained after introducing the ventilation coefficients f v for water vapor diffusion and fh 
for heat transport in Equation (3.9). That is 


r = (l + 


fvsj^S 

KC p T 2 


pgy. 

\np l A i f v N p r v 


(3.12) 


where N p and r p are the concentration and the radius of precipitating ice particles, 
respectively. Let r p =10 3 pm, A p =10 3 m°, T = 260 K and p = 550 hpa. The timescale of 
water vapor deposition on snowflakes is 1184 s when f v = 2, and the timescale of 
deposition on graupel is 236 s when/ v =10. Since both timescales are much larger than 10 
s, water vapor deposition on to snowflakes and graupel is suitable for explicit simulation 
in cloud-resolving modeling. 


c. Ice fusion 

An air parcel with uniform spherical ice particles is studied for the timescale of cloud 
ice fusion. For simplicity, it is assumed that no liquid water surrounds ice particles. Let A, 
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and Yi denote the concentration and the radius of ice particles, respectively. The melting 
of an ice particle is described with (e.g., Pruppacher and Klett 1997). 


- AnpL f r^ ^ = 4 mf h K a (T - T) + AnrJL (— — E ™ {T o l 
y, f , d{ ,jh ot ,Jy ^ R T R Jo 


), 


(3.13) 


where T 0 =273.15 K and e is the partial pressure of water vapor in the air. In the preceding 
equation, the term on the left side represents the latent cooling due to ice fusion, and the 
two terms on the right side represent the sensible and latent heat fluxes from air to the ice 
particle, respectively. Obviously the latent heat of water vapor condensation at the 
surface is balanced by part of the latent heat of ice melting, and the latent heats do not 
change air temperature directly. 

Assume that the parcel is adiabatic and stationary. For ice fusion, heat is transfered 
from air to ice particles. Thus air temperature is decreased, which is described as 


C,^- = -4*rthK.^V-T.)- 

dt p 


The contrast between Equations (3.14) and (2.1) shows the timescale of ice fusion 

pc. 


(3.14) 


T = 


4 *f h K a N iTi 


(3.15) 


Using the above expression, the timescales of ice fusion and snowflake fusion are 
estimated with fh=\ and 2, respectively. Their results are displayed in Figure 4 with thin 
and thick lines, respectively. 

Equation (3.15) and Figure 4 are also suitable for the freezing of water that is 
collected by ice particles. As shown in Figure 4, the timescales in relation to fusion and 
freezing are larger than 10 s. Therefore, the fusion and freezing processes are suitable for 
explicit simulation in cloud-resolving modeling. 
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4. Moist entropy as a prognostic variable 

The microphysical timescales analyzed in the preceding section show that only the 
timescale of water vapor condensation is smaller than the typical time step used in cloud- 
resolving modeling. Since this short timescale may bring about computational 
phenomena in cloud-resolving modeling, the choice of prognostic thermodynamic 
variables is discussed in this section for the removal of the computational phenomena. 

a. Prognostic thermodynamic variables 

To remove the computational phenomenon in Equation (2.2) when A t>t , it is 
assumed that 0=0 in one step. As shown in Figure 1, 0=0 approximates the analytical 
solution well when A t>r. In other words, the prognostic variable 0 degenerates into a 
diagnostic variable. Next, prognostic thermodynamic variables are discussed for cloud- 
resolving modeling. 

In current cloud-resolving models, three prognostic thermodynamic variables are used 
to simulate non-precipitating warm clouds and those three variables plus others to 
simulate other clouds (e.g., Grabowski 1989, Tao and Simpson 1993, Tompkins and 
Craig 1998, Tao et al. 2003b). The three prognostic variables are T temperature (or its 
equivalent), q v the mixing ratio of water vapor and q c the mixing ratio of cloud water. 
However the three prognostic variables, as shown in Equations (3.4) and (3.5), involve 
one timescale that is smaller than 10 s the order of time steps used in cloud-resolving 
modeling. Since only one short timescale is involved, it is impossible for all three 
prognostic variables to degenerate into diagnostic variables at the same time. Thus, it is 
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suggested that other prognostic variables (e.g., moist entropy) be used for cloud-resolving 
modeling. 

Moist entropy per unit mass of dry air is expressed as (Zeng, Tao and Simpson 2004) 
s = ( C p + c i<h ) In ~ R d In — ~ + ^q v -R v q v \nf-^ r q i , (4.1) 

T ref Pref T T 

where relative humidity f=e/E sw \ T re f = 273.15 K and p re f = 10 5 pa are the reference 
temperature and pressure, respectively; and the total mixing ratio of airborne water (water 
vapor, cloud water and ice) 

<7, = <7v + <7c + <7, • ( 4 - 2 ) 

For the simulation of non-precipitating warm clouds, the three prognostic variables T, 
q v and q c can be replaced equivalently with three other prognostic variables s, q, and the 
supersaturation ( qjq vsw -1) . Of all the other three variables, as shown in Equation (3.6), 

only the supersaturation involves the short timescale of water vapor condensation. Since 
cloud condensation (or evaporation) can not change 5 or q t , it is reasonable to assume that 
the supersaturation is zero when cloud water exists 1 . Consequently, three prognostic 
variables are decreased to two prognostic variables (i.e., 5 and q,). 

Zero supersaturation implies that air is saturated with respect to water and real 
supersaturation is removed. In fact, it is impossible for real supersaturation to be 
simulated in cloud-resolving models even though a small time step is used, because no 
spectra of cloud particles are represented explicitly. 

1 In spectral-bin models (e.g., Tao et al. 2003b), the spectra of cloud particles are represented 
explicitly. Thus, the approximation of zero super saturation can be removed and real 
supersaturation can be simulated properly with a small time step. 
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Fortunately, the approximation of zero supersaturation is supported by observations. 
Since real supersaturation in clouds can not be measured directly, it is usually calculated 
with observational data of temperature, vertical velocity and cloud droplets. It is found 
that real supersaturation in cumulus clouds is in the range from -0.5 to 0.5% and rarely 
exceeds 1% (Politovich and Cooper 1988). 

The approximation of zero supersaturation is understandable in the atmosphere. When 
no ice is involved, air containing cloud water is very close to saturation over water 
because the concentration of cloud condensation nuclei is very high (e.g., Howell 1949; 
Politovich and Cooper 1988). When ice particles exist, air with cloud water is still close 
to saturation over water because of the great difference of the concentrations of cloud 
drops and ice particles (Korolev and Mazin 2003, Zeng, Tao and Simpson 2004). The 
concentration of cloud drops is about four orders larger than the concentration of ice 
particles (e.g., Wallace and Hobbs 1977, Pruppacher and Klett 1997). Thus the timescale 
of cloud water evaporation is much smaller than the timescale of water vapor deposition 
on ice particles (see Figures 2 and 3 for comparison). As a result, the air with cloud water 
is almost saturated with respect to water, which provides an environment for the 
Bergeron process in cold clouds (Bergeron 1935). Summarily, it is physically clear that 
( s , q t ) can replace (T, q v , q c ) as prognostic variables in cloud-resolving modeling. 

b. Diagnosing temperature from moist entropy 

This subsection describes how to diagnose temperature from moist entropy and other 
prognostic variables. The procedure is summarized in the flow chart in Figure 5. When p 
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the pressure and q t the mixing ratio of cloud ice are known, ( T , q v , q c ) the air temperature 
and the mixing ratios of water vapor and cloud water can be diagnosed from (s, q,') the 
moist entropy and the total mixing ratio of airborne water. 

Consider an imaginary parcel with s the moist entropy , p the pressure and q, the 
mixing ratio of cloud ice. The imaginary parcel stays right at water-saturation and 
contains no cloud water. Its temperature is denoted as T w * and its saturation mixing ratio 
of water vapor as 


E (T > 

= 0.622- w V wJ 


E (T ) 

sw V w J 


(4.3) 


Using the preceding equation. Equation (4.1) is solved for T w * with the Newton iterative 
method first. Then, q V5W * is determined by Equation (4.3). 

Since moist entropy and the total mixing ratio of airborne water are conserved when 
water vapor condensates or cloud drops evaporate, q vsw for the imaginary parcel is 
compared with q r qt for the original parcel, showing whether the original parcel is 
saturated with respect to water. When q,-qi<q V sw*, the air is unsaturated with respect to 
water. Thus q c = 0, q v =q,-q„ and Equation (4.1) is solved for the air temperature T with the 


2 When cloud ice is involved, moist entropy changes slightly with temperature in adiabatic 
motion. Thus, the moist entropy of the imaginary parcel is adjusted to 

s + q, [L f (T)/T — Lj (T*)/T* ] so that the imaginary parcel can change adiabatically to the state 
where moist entropy is the same as that of the original parcel [see Equation (3.5) of Zeng et al. 
(2004) for details]. In other words, s + q t L f (T)/T is conserved in cloud condensation (or 
evaporation). 
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Newton iterative method. When q r q,>q vsw *, the air is saturated with respect to water. As 
a result, q v — qvsw i T—T w and q ( — q i~q v~q i - 


c. Comparison of two numerical models 

In this subsection, two numerical models are used to simulate an air parcel in adiabatic 
upward motion. They use ( T , q v , q ( ) and (s, q t ) as prognostic variables, respectively. Their 
results are compared to show the benefit of using moist entropy as a prognostic variable. 

In the first model, (T, q v , q c ) are used as prognostic variables. Their governing 
equations are 


dq c _ dq v 

dt 


dt 


.n . \ dT „.^ dln(p-g) _ d lne 

(c p +c pv q v c,q c ) dt i dt q v J dt 


dp 

- = -pgw 

dt 


= -L 


dq v 

v dt 


(4.4) 

(4.5) 

(4.6) 


and 


dfL 

dt 


4vsw.L v _\- 1 Qv 

2 ' 


^ *3 VSW^V ^ -1 *3 y *3 V5W 


KC P T 


(4.7) 


when q v >q vsw - When q v <q vsw , Equation (4.7) still works if q c > 0. Otherwise, dq v ldt= 0. 

Of the preceding governing equations, Equation (4.6) is obtained from the hydrostatic 
equation, and Equation (4.7) is obtained from Equations (3.4) and (3.8). As shown in 
Equation (3.8), the timescale t changes with time. Its variation can be simulated 
explicitly in spectral-bin models that explicitly represent the spectrum of cloud drops 
(e.g., Tao et al. 2003b). For simplicity, a constant timescale T=l s is used in Equation 
(4.7) so as to bring the computational phenomena into focus. 
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In the second model, (s, q,) are used as prognostic variables. Their governing 
equations are dst dt-dq x ! dt - 0 and Equation (4.6). Both models are equivalent. They 
involve neither precipitation nor ice. They follow the numerical scheme in Equation 
( 2 . 2 ). 

The air parcel is assumed to move upward with a vertical velocity of w = 4 m/s as well 
as an initial pressure 1000 hpa, relative humidity of 85% and temperature of 300 K. The 
numerical results from the two models are displayed in Figure 6. The results from the 
first model with Ar=0.1 s are displayed with thin solid lines. These results can be 
regarded as a benchmark to check the results in other cases. 

The results from the first model using Ar=3 s are displayed with thin dashed lines in 
Figure 6. As shown in the figure, the results are bad compared to those with Ar=0.1 s, and 
both spurious supersaturation and negative mixing ratio of cloud water are present. When 
the model uses Ar=10 s, it blows up due to computational instability. The results from the 
second model with At=10 s are displayed with thick dashed lines in Figure 6. As shown 
in the figure, the results agree well with those of the first model with A?=0.1 s, and 
neither spurious supersaturation nor negative mixing ratio of cloud water appears. These 
numerical tests clearly show that moist entropy and the total mixing ratio of airborne 
water are more efficient than temperature and the mixing ratios of water vapor and cloud 
water as prognostic variables in cloud-resolving modeling. 

5. Summary 

Spurious computational phenomena occur in cloud-resolving models when the time 
step for explicit integration is larger than a microphysical timescale. The phenomena are 
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evident as spurious supersaturation and negative mixing ratio of cloud water (e.g., 
Grabowski 1989, Tao et al. 1989). The related computational error can affect long-term 
cloud-resolving modeling significantly (e.g., Robe and Emanuel 1996). 

In this paper, the microphysical timescales in clouds are studied in comparison with 
the time step for cloud-resolving modeling, providing theoretical evidence for the choice 
of prognostic thermodynamic variables. Of all the microphysical timescales that are 
analyzed in the paper, only the timescale of water vapor condensation (or cloud water 
evaporation) is smaller than 10 s, the order of a typical time step used in cloud-resolving 
modeling, indicating that cloud condensation (or evaporation) is not suitable for explicit 
simulation in cloud-resolving modeling. 

In current cloud-resolving models, the three variables (T, q v , q c ) are used as prognostic 
variables to simulate non-precipitating warm clouds. Since all three variables involve the 
short timescale of water vapor condensation, it is suggested that they be replaced with 
moist entropy s, the total mixing ratio of airborne water q, and supersaturation. The 
supersaturation, rather than 5 and q t , involves the short timescale. Therefore, it is 
reasonable to assume zero supersaturation when cloud water exists, and consequently the 
three prognostic variables are decreased to two prognostic variables (i.e., 5 and q t ). 

Two numerical models are used to simulate an air parcel in adiabatic upward motion. 
Their results are compared to show the benefit of using moist entropy as a prognostic 
variable. The first numerical model uses ( T , q v , q c ) as prognostic variables. When its time 
step is small (e.g., less than 1 s), its results are reasonable. When its time step is large 
(e.g., larger than 2 s), its results are overwhelmed with great computational errors, 
exhibiting the phenomena of spurious supers aturati on and negative mixing ratio of cloud 


17 



water. The second model uses (s, q t ) as prognostic variables. When its time step is 10 s 
(or even larger), its results still agree well with those of the first model using a very small 
time step (e.g., 0.1 s). The comparison of the two models shows clearly that (s, q t ) can 
work efficiently as prognostic variables in cloud-resolving modeling. 
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APPENDIX A 


List of Symbols 

Aw/Ai : defined in (3.3)/(3.10) 

Cp/Cpv/ci : specific heat of dry air/water vapor/liquid water 

D v : coefficient of water vapor diffusion in air 

e : partial pressure of water vapor 

Esw/Esi : saturation vapor pressure over water/ice 

f=e/E sw : relative humidity 

f/fh : ventilation coefficient for water vapor diffusion/heat transport 

g : acceleration due to gravity 

K a : coefficient of air heat conductivity 

Ly/L/Lf : latent heat of vaporization/sublimation/freezing 

NJNr/Ni/Np : concentration of cloud droplets/raindrops/cloud ice particles/precipitating 
ice particles 

p : total pressure of moist air 
Pre/= 10 5 pa : reference pressure 

qvsdqvsi '■ saturation mixing ratio of water vapor over water/ice 

qv/q/q/q/q/qg : mixing ratio of water vapor/cloud water/cloud ice/rain/snow/graupel 

q t =q v +q c +qi : total mixing ratio of airborne water 

r c lr r lrilr p : radius of cloud droplets/raindrops/cloud ice particles/precipitating ice particles 
R/R v : gas constant of dry air/water vapor 
s : moist entropy per unit mass 
t : time 
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T : temperature 

r o =273.15 K : absolute temperature at the melting point 
T re f =273.15 K : reference temperature 
w : vertical velocity 
z : height 
r: timescale 

At : time step for integration 
p : air density 

pjpi : density of liquid water/ice 
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Figure Captions 


Figure 1 The variable <j> in the vertical axis varies with time when At=10 s and the 
timescale ris 2, 5 and 12.5 s, respectively. The numerical and the analytical solutions are 
shown with thick and thin lines, respectively. 

Figure 2 The timescale of water vapor condensation varies with pressure (or 
temperature). The results with N c r c = 500 and 3000 pm-cm' 3 are displayed with thin and 
thick lines, respectively. 


Figure 3 The timescale of water vapor deposition varies with pressure (or temperature). 
The results with A^rplO' 3 , 10' 1 and 10 1 pm-cm' 3 are displayed with dashed, thin and thick 
lines, respectively. 


Figure 4 The timescale of ice fusion varies with N p r p . The timescales for cloud ice 
fusion and snowflake fusion are displayed with thin and thick lines, respectively. 


Figure 5 Schematic on the diagnosis of temperature from moist entropy and the total 
mixing ratio of airborne water. 


Figure 6 Change in output variables with height for two different models. The variables 
in a traditional model with Ar=0.1 s are displayed with thin solid lines; those in the 
traditional (or first) model with At=3 s are displayed with thin dashed lines; and those in 
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the present (or second) model with Ai=10 s are displayed with thick dashed lines. Thin 
solid lines coincide with thick dashed lines. 
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Figure 1 The variable (p in the vertical axis varies with time when Ar=10 s and the 
timescale ris 2, 5 and 12.5 s, respectively. The numerical and the analytical solutions are 
shown with thick and thin lines, respectively. 
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Figure 2 The timescale of water vapor condensation varies with pressure (or 
temperature). The results with N c r c = 500 and 3000 pm-cm' 3 are displayed with thin and 
thick lines, respectively. 
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Figure 5 Schematic on the diagnosis of temperature from moist entropy and the total 
mixing ratio of airborne water. 
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Figure 6 Change in output variables with height for two different models. The variables 
in a traditional model with Ar=0.1 s are displayed with thin solid lines; those in the 
traditional (or first) model with Af=3 s are displayed with thin dashed lines; and those in 
the present (or second) model with Ar=10 s are displayed with thick dashed lines. Thin 
solid lines coincide with thick dashed lines. 
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Popular Summary 

Spurious computational phenomena (e.g., spurious supersaturation and negative mixing 
ratio of cloud water) exist in current cloud-resolving modeling. They originate in the 
interaction between model dynamics and cloud microphysics. In this paper, theoretical 
analysis is given to show that the phenomena appear while the timescale of cloud 
condensation or evaporation is smaller than the time step used in cloud-resolving models. 

To remove the phenomena, it is suggested that moist entropy be used as a prognostic 
variable and temperature be diagnosed from that and other prognostic variables. Numerical 
tests show that the spurious computational phenomena, as expected, disappear when moist 
entropy is used as a prognostic variable. 


